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ABSTRACT 

We identify close companions of Brightest Cluster Galaxies (BCGs) for the pur- 
pose of quantifying the rate at which these galaxies grow via mergers. By exploiting 
deep photometric data from the CFHTLS, we probe the number of companions per 
BCG (N c ) with luminosity ratios down to those corresponding to potential minor 
mergers of 20:1. We also measure the average luminosity in companions per galaxy 
(L c ). We find that N c and L c rise steeply with luminosity ratio for both the BCGs, and 
a control sample of other bright, red, cluster galaxies. The trend for BCGs rises more 
steeply, resulting in a larger number of close companions. For companions within 
50kpc of a BCG, N c = 1.38±0.14 and L c =2.14±0.31x 10 10 L o and for companions 
within 50kpc of a luminosity matched control sample of non-BCGs, N c =0.87±0.08 
and L c =1.48±O.2OxlO lo L . This suggests that the BCGs are likely to undergo more 
mergers compared to otherwise comparable luminous galaxies. Additionally, compared 
to a local sample of luminous red galaxies, the more distant sample presented in this 
study (with redshifts between 0.15-0.39,) shows a higher N c , suggesting the younger 
and smaller BCGs are still undergoing hierarchical formation. Using the Millennium 
Simulations we model and estimate the level of contamination due to unrelated cluster 
galaxies. The contamination by interloping galaxies is 50% within projected separa- 
tions of 50kpc, but within 30kpc, 60% of identified companions are real physical 
companions. We conclude that the luminosity of bound merger candidates down to 
luminosity ratios of 20:1 could be adding as much as 10% to the mass of a typical 
BCG over 0.5 Gyr at redshifts of z~0.3. 
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1 INTRODUCTION 

Deducing the formation and evolution of galaxies is an active 
area in both observational and theoretical work in astron- 
omy. Hierarchical formation of galaxies through mergers is 
certainly among the main processes in a galaxy's evolution- 
ary history. In modern studies, close galaxy pairs are used as 
a diagnostic of merger activity, as they should be precursors 
to mergers (e.g., |Kitzbichler &; White 2008). 

Galaxies with close companions have been found to have 
higher asy mmetry (|Patton et alT||2000| | Hernandez- Toledo" 



et al.|2005| [Patton et al.|2005||De Propris et al.|2007||McIn- 

tosh et al. 2008|, e nhanced star formation (|Lambas et al 
20031 lAlonso et al. 



2007 



lEllison et al. 



20041 [Nlkolic et al.||2004( |Barton et al. 
2008) |Li et al.|2008| [ Ellison et aI.|20ToT 



Woods et al. 2010 ), bluer central colours (|Patton et al 
lower metallicities (Kewley et al. 2006; Ellison et al 



2011), 



2008) 



and higher AGN fractions ( Ellison et al.||20"TT ) than galax- 



ies without close companions. This indicates that mergers 
can significantly affect the morphological, star forming and 
chemical properties of the galaxies involved in these inter- 
actions. 

The frequency of close galaxy pairs can be used to deter- 
mine the merger fraction, and this has been calculated from 
numerous studies using large galaxy surveys at various red- 



shifts (e.g.,|Patton et al.|2002||Kartaltepe et al.|2007 


Hsieh 


et al.|2008; Lin et al.||2008| |Patton & Atfleld| 2008 ; 


Rawat 


et al. 2008 


|de Ravel et al.|2009| |Lopez-Sanjuan et al.|2010| 


Lotz et al. 


201 1|. The merger rate depends on galaxy mass, 



the luminosity ratio between merging galaxies and the sur- 



rounding environment \ Khochfar &; Bur kert 200 1 | |Berrier 
eTaLl|2006l |Cox et al.||2008| |Guo fe White||2008| |Patton 
Atfie ld|2008 ). No one environment is host to merging galax- 
ies, as galaxy groups ( Goto|2005 Brough et al.|20 06; Nolan 
et al.||2007 ), cluster infall regions ( van Dokkum et al.| |1999 
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Tran et al.|200 5; Moss 2006), and field galaxies (Barton et al. 
2007) have all been shown to undergo mergers. Ellison et al. 



results from dark matter simulations of IDe Lu cia & Blaizot 



( 2010 ) find that lower density environments contain a higher 
fraction of pairs, but high density environments such as clus- 



ters still contain a significant population of pairs ( |Masjedi 
et al.|2006| jPatton et al.|2011| ) . What is more, major merg- 



ers in luminous cluster galaxies continue to occur, even at 
low redshift ( |McIntosh et al.|2008| [Brough et aL|2011 ). 

The period of time between z=0-l has seen a growth 



of the red sequence (Wake et al. 2006 |Faber et al. 2007 



Brown et al. 2007 Ruhland et al. 2009) which is thought 



to be spurred by major merging events ( Schiminovich et al. 
2007 Wild et al.|20 09), at least in L* galaxies. However, the 



importance of minor mergers on the build up of the red se- 



quence has recently gained attention. Lotz et al. (2011) find 



that the rate of minor mergers (10:1) is three times that 



for major mergers (4:1). Shapiro et al. ( 2010| have studied 
the morphology and kinematics of a large sample of ellip- 
tical type galaxies and conclude that recent star forming 
events are likely to be the product of minor mergers that 
are building up the galaxy's bulge. Similarly, Kaviraj et al.| 



(2010) find too many morphologically disturbed and star 



forming ellipticals to be accounted for by major mergers 
alone. Since minor mergers are predicted to be a more com- 
mon phenomena ( Guo &; White|2008 ), it follows that minor 
mergers largely influence the build up of elliptical galaxies 
at late epochs. 

In the hierarchical view of galaxy formation, the most 
massive galaxies would be those most recently formed. At 
first, this seems at odds with the idea of downsizing, where 
the most massive galaxies formed first. Understanding the 
formation and evolution of these most massive galaxies, 
therefore, is of particular importance. The most massive 
galaxies are Brightest Cluster Galaxies (BCGs), where the 
bulk of the stellar mass is from extremely red and very old 
stars. BCGs live in a particularly complex environment at 
the gravitational potential centre of clusters, where dense X- 
ray emitting intracluster gas is often found. The galaxies in 
these massive X-ray clusters may be old simply because they 
were the first density peaks of the early universe to collapse 



(Treu et al. 2005). However, their formation and evolution 



might not be straightforward. Particularly intriguing are in- 
dications of recent star formation, such as blue colours and 
optical line emission which have been observed in BCGs that 
lie at the X-ray peak of cool -core clusters (|Crawford et al 



1999; Mc Namara et al. || 2 004 1 [Edwards et al.||2007| |Rafferty 



et al.||2008| jBildfell et al. 2008). Thus, cluster-scale prop- 
erties likely factor into the presence of emission lines seen 
in BCGs, but other triggers may exist. Outside of cool-core 
clusters and BCGs, 10% of luminous red galaxies show op- 
tical emission lines at their centres, regardless of whether 
the local galaxy density describes a massive cluster, or poor 
group ( |Edwards et al.||2007| ). Could this conspicuous line 
emission in BCGs and other luminous red galaxies be the 
result of recent mergers? The cool gas associated with cool- 
core BCGs may be influenced by interactions. |Donahue et al.| 
(2007} and |Wilman et al.] ( |2006| point out that galaxies with 
more disturbed Ha morphologies have companions and ex- 
plicitly ask whether galaxy-galaxy interactions are required 
to trigger the central nebulae found in cooling core BCGs. 

Evidence that recent and current galaxy interactions 
are an important component of BCG evolution includes the 



(2007) which show that half of the BCG assembly has taken 
place since z^ 0.5. Observational evidence for this is also 
beginning to be compiled: Galaxy kinematics from [Brough 
et al. (2011) have recently confirmed that nearly equal mass 



pairs of two different local BCGs are gravitationally bound 
and likely to merge within 0.3 Gyr. More generally, in a sam- 
ple of 515 BCGs, [Liu et al.| ( [2009| find that dry major merg- 
ers are an important component of the buildup from z=0.12 
to the present, but that not all the mass of the secondary 
is incorporated into the merger. Much of recent growth in 
massive galaxies may be in lower mass-ratio mergers. For 
example, |Stott et al.| ( |2011| find that BCGs have increased 
in size by only 30% between z=0.25 and 1, and although 



Bernardi (2009) and Ascaso et al. (2011) find a significant 
increase in BCG size between z=0.25 and 0, the BCG lumi- 
nosity profile has not changed, suggesting gentle processes 
like minor mergers. In the theoretical work of |Naab et al.| 
(2009) minor mergers are also found to be an important 
component of the size evolution of massive galaxies, and re- 
cent observations of distant galaxies by Bluck et al. (2012) 



and more local systems |Trujillo et al. (2011) support these 
findings. 

In this paper, we focus on close companions to BCGs in 
massive galaxy clusters. The large mass and high luminos- 
ity of BCGs allow for large luminosity ratios between the 
parent galaxy and companion. We use the deep photome- 
try of the Canada- France- Hawaii Telescope Legacy Survey 
(CFHTLS) program to study the galaxies, so that we can 
probe BCGs with 20:1 luminosity ratio companions; poten- 
tially major and minor mergers. Specifically, we measure 
the number (N c ) and luminosity (L c ) of the companions 
and their proximity. Using a comparison sample of lumi- 
nous galaxies gathered from the Millennium simulations, we 
estimate the level of contamination and calculate the dy- 
namical friction time-scale for close companions near BCGs. 
The paper is organized as follows: in Section [2] we discuss 
our datasets and method for calculating N c and L c , in Sec- 
tion [3] we present our results and how they may be affected 
by interloping galaxies, in Section 4 we discuss the relative 
importance of minor and major mergers on BCG and non- 
BCG evolution, and in Section[5]we present our conclusions. 
Throughout the paper we adopt a cosmology with Q m =0.3, 
Q A =0.7 and H Q =71 krns" 1 Mpc -1 . 



2 DATA AND METHOD 

2.1 CFHTLS BCGs and controls 

The BCGs and control galaxies are chosen from a list of 
clusters found in the CFHTLS wide fields (Wl, W2, W3 and 
W4) data release version T0004. This release is complete for 
g', r', and 1 in target depth, though it covers only 119 of 
170 square degrees on the sky. The survey depths are 25.5, 
24.8, and 24.5, in AB magnitudes for the 3 colour bands, 
respectively. We use the clusters from|Lu et al.|(|2009| f ound 
using the red sequence method of |G ladder s & Yee ( 2005 ) . All 
redshifts in this study are the photometric redshifts derived 
from their analysis. There are 127 clusters in their study 
for which there are at least 12 red members brighter than 
M*+2 magnitudes, within 0.5 Mpc of the cluster centre, and 
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Figure 1. Top - For illustrative purposes, we show five- colour images of four representative BCGs in our sample of 91. They display 
relaxed, regular morphologies, typical of BCGs. From left to right, the BCGs have (M r ,z) of (-22.86,0.25), (-22.96,0.25), (-23.45,0.24), 
(-23.59,0.28). Bottom - Four galaxies from the sample of 197 luminosity-matched controls. From left to right, the BCGs have (M r ,z) of 
(-22.85,0.23), (-23.01,0.25), (-23.34,0.23) ,(-23.18,0.28) . All figures are 120" on a side, which corresponds to 425kpc at z=0.25, where 
z is the photometric cluster redshift derived from the red sequence fit. 



which have at least 50 galaxies within 1 Mpc of the cluster 
centre. Figures 7 and 9 of Lu et al. ( 2009| show stacked 
colour magnitude diagrams from the cluster sample, as well 
as the colour magnitude diagram for the richest cluster in 
the sample, respectively. 

We search for potential cluster members within 1 Mpc 
of each cluster centre. We use the publicly available pho- 
tometric catalogues (from TERAPIX, and available on- 
line through the Canadian Astronomical Data Centre) and 
gather all available optical bands from the full T0006 data 
release, the first complete version of the catalogues including 
u* and z' photometry (with respective AB depth of 25.3 and 
23.6). Because the clusters cover a range of redshifts between 
z=0. 15-0.39, we apply the K-correction code of Blant on~fc] 
Roweisl (12007). 



Lu et al. 



(2009), we 



Pipino et al. (2011) 



We pick the BCG as the brightest galaxy in r' which is 
red ((u*-r / )>2.22) and also within 0.25 Mpc of the cluster 
centre, there are 101 of these. Because this is essentially 
based off of the red-sequence catalog of 
are only looking at red BCGs. Recently 
has found that blue BCGs make up a small fraction of the 
BCG population for clusters at redshifts between 0.1 and 
0.3; our sample does not contain these blue BCGs. 

We construct a volume limited survey of BCGs with 
absolute magnitudes brighter than M r =-20.06 (r'=19.2 at 
z=0.15). By nature, BCGs are at the bright end of the galaxy 
luminosity function, with typical magnitudes of M r =-23, 
thus our magnitude cuts well sample the BCG population. 
As the faintest BCG in the sample is actually brighter than 
our cut by 1.54 magnitudes, our sample essentially includes 
companions down to magnitude differences of 4.8. However, 
we choose to restrict our analysis to companions with mag- 
nitude differences of up to 3.25 (20:1 companions) since the 



purity of the sample degrades when looking at fainter com- 
panions, as we will see in Section |cT4| 

The sample of BCGs was examined by eye using the 
CFHTLS cutout service ( |Gwyn|2008"l ). We found that many 
of the bright sources were in fact stars that had not been 
flagged. We re- flagged these in our galaxy tables and re- 
ran the BCG finder, there was one case where a star was 
too close to the potential BCG to accurately measure its 
magnitude, so we deleted this galaxy from our sample. It 
can happen that an obvious foreground elliptical can satisfy 
these criteria, and in our by-eye inspection we found 9 such 
cases. The final sample contains 91 BCGs. The five-colour 
postage stamps of four of the BCGs taken from the cutout 
service are shown in Figure [I] Figure [2] shows the distribu- 
tion of cluster redshifts, (u*-r ; ) colour, and M r for the BCG 
sample. The BCG population has a distribution of redshifts 
which peaks at z^0.3, (u*-r')~3.2, and bright magnitudes 
of M r — 22.5. 

We construct a luminosity matched control sample 
which we use to compare to the cluster BCGs. We accom- 
plish this by collecting up to 3 other luminous red galaxies 
(LRGs) near the cluster centre, excluding the BCG, and in- 
cluding only galaxies with luminosities similar to the BCG. 
We also include only red galaxies on the same red sequence 
as the host BCG by using colour cuts of (g'-r^zb 0.075 and 
(r'-i^iO^ ( Lu et al.|2009 ). Again, all controls are inspected 
by eye and any lingering bright stars are removed. Figure [I] 
shows four of our 197 controls. Figure [2] shows the redshift, 
colour, and absolute magnitude distributions of the control 
sample. It is well matched to the BCG sample in redshift 
and absolute magnitude. The colours of the controls are not 
as heavily weighted towards the red as in the BCGs because 
the latter were chosen to be the reddest. 
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Figure 2. Parent galaxy properties: The distribution of galaxy properties for the sample of 91 BCGs (solid line) and 197 controls 
(dashed line) are shown. Left - The distribution of redshifts; the full redshift range has been sampled, centre - The distribution of parent 
galaxy colour. The BCGs are chosen to be the reddest and most luminous of the cluster, and the sample is skewed towards redder colours. 
The control sample is more evenly distributed. Right - The distribution of absolute M r magnitudes assuming the cluster redshift. By 
definition these are massive galaxies. Both distributions peak at M r =-22.5. 



2.2 Identifying close companions 

Companions to each BCG and control galaxy are recovered 
by searching all galaxies in the photometric tables (red and 
blue) within 30kpc and 50kpc of the BCG, using the red- 
shift derived by |Lu et aL ( 2009 ) . The strongest evidence 
for interactions and mergers in galaxy pairs has been found 
within 30kpc and this smaller search radius results in sys- 
tems which are more likely to merge ( |Patton et al.|[2011| 
Wong et al. |2011| ). However, the larger search radius natu- 
rally leads to a larger sample of companions. Additionally, 
because the BCGs are the most massive systems, the larger 
radius of 50 kpc may be more appropriate. 

We calculate the galaxy luminosity assuming the clus- 
ter redshift derived from the red sequence fitting. We have 
verified that none of the close companions to the BCGs is 
a control, and likewise, none of the control galaxy compan- 
ions is a BCG. The combination of selecting the BCG sam- 
ple, which is intrinsically bright, with the deep photometry 
obtained with the CFHTLS therefore allows us to explore 
properties of close companions, ratios of 2:1, 4:1, 10:1, and 
down to 20:1, larger magnitude ratios than previously ex- 
plored. These galaxies are those which will potentially lead 
to both major and minor mergers. In Section|3]4]we show the 
relationship between the companions we identify based on 
the projected separation to the parent galaxy, versus those 
which are close in three dimensions. Section [2.3| which fol- 
lows, introduces the simulation data we use to find the 3D 
separations. 



2.3 BCGs in the Millennium Simulations 

We use output from the Millennium Simulations, where 3D 
separations can be measured, in order to quantify the pu- 



rity of the CFHTLS results. We use the |Blaizot et al.| (|2005) 
mock redshift catalogue to do this. This catalogue contains 
over 5.7 million galaxies, each with apparent magnitudes, 
calculated stellar and virial masses, redshifts, and positions. 
To find a sample of simulated galaxies which mimics lu- 
minous red galaxies, we choose the ^300000 with M r <-24. 



However, most of these luminous galaxies are not at the 
gravitational potential of a massive cluster. Therefore, to 
better match to the observed data from the CFHTLS, of 
BCGs and controls which are embedded within massive clus- 
ters, we construct a second sample of potential CFHTLS 
BCG-like galaxies, and then collect the companions in the 
same way as for the observed data. There are 291384 sim- 
ulated red (u / -r / >2.22) and luminous (-24.5<M r <-21.8) 
galaxies within the redshift range 0.2-0.4. From this list, 
we keep only those that are the brightest within a radius of 
0.5 Mpc, resulting in 272042 galaxies. These bright galaxies 
have a total of 138395 companions, within 0.5 Mpc. 

As alluded to by the previous numbers, the vast major- 
ity of these galaxies, especially those with the lower intrin- 
sic luminosities, are likely to be luminous galaxies in poor 
groups, or at the edge of clusters. To select BCGs that most 
accurately reflect those of the CFHTLS sample, we count 
the number of luminous red galaxies within 2 magnitudes 
of the potential BCG and within 0.5 Mpc. Only if there are 
at least 3 other luminous red galaxies with velocities within 
1500 km/s of the BCG, we deem the potential BCG to be 
part of a massive cluster, 260 of these CFHTLS BCG-like 
galaxies remain. To find the companions to these BCGs, we 
search the full |Blaizot et aL] ( |2005| catalog for all galaxies, 
red and blue, within 0.5 Mpc of the BCG. There are 871 
companions. 



3 RESULTS 

To measure the importance of companions for BCGs and 



the controls, we follow the methods of Patton et al. (2000). 
For companions within a projected distance of 30 kpc and 
50 kpc of the host galaxy, we count the number of galaxies, or 
add up their luminosities, and divide by the total number 
of parent galaxies, yielding the average number (N c ) and 
luminosity (L c ) of companions. 

There are 126 BCG companions with luminosity ra- 
tios down to 20:1 and within 50 kpc (46 within 30 kpc,) 
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Figure 3. Left - N c as a function of the luminosity ratio for 
companions to BCGs (black) and controls (grey). Right - L c as a 
function of the luminosity ratio for companions to BCGs (black) 
and controls (grey). Companions within SOkpc are shown on the 
top. Those within 50 kpc are shown on the bottom where BCGs 
have higher values of N c and L c than the control sample. The last 
point includes the impure set of all measured companions beyond 
magnitude ratios of 20:1. 



and 172 control companions (77 within 30 kpc). All of 
the companions have been inspected by eye, and obvious 
stars removed. The average N c for BCGs is 1.38±0.14, 
larger than for the luminosity-matched control sample, 
which has N c =0.87±0.08. The average L c of BCGs, 
2.14±0.31xl0 10 L©, is also higher than that for the non- 
BCG sample, 1.48±0.20x 1O 1O L , where the uncertainties 
are calculated from the standard error on the mean. At face 
value, these results indicate that BCGs accrete more com- 
panions and more stellar mass more often that non-BCGs. 
We will return to this issue after addressing contamination 



by interlopers (no n- merging companions) in Section 3.4 



3.1 The number and luminosity of companions as 
a function of luminosity ratio 

We construct the cumulative distribution of N c and L c as a 
function of the luminosity ratio between the parent (BCG 
or control) and the companion. 

Figure [3] shows that the BCGs have similar low numbers 
of 2:1 and 4:1 companions as do the controls. The dearth of 
near mass companions is expected as these luminous clus- 
ter galaxies are at the bright end of the cluster luminosity 
function and rare, thus less likely to have a close companion 
with a similar size. The trends for BCGs and controls are the 
same to within errors for the sample of companions within 
30 kpc. It is interesting, however, that the number of com- 
panions for BCGs rises more steeply than for the controls, 
and is consistently higher. In the sample that is statistically 



Figure 4. Left - N c as a function of the magnitude of BCGs 
(black) and controls (grey). Right - L c as a function of the mag- 
nitude of BCGs (black) and controls (grey). The sample with 
r p <30 kpc is shown on the top, and the sample with r p <50 kpc 
is shown on the bottom. All panels include companions down to 
luminosity ratios of 20:1. 



more significant, including companions out to 50 kpc, the 
number of 10:1 and 20:1 companions is larger for the BCGs 
than for the controls. 

This steeper rise for companions of BCGs can also be 
easily seen for L c (Figure [3] right panels). There is more 
luminosity in companions per BCG than in companions per 
control when integrating over all companions found, includ- 
ing those with ratios beyond 20:1. 



3.2 The number and luminosity of companions as 
a function of parent galaxy absolute 
magnitude 

The number of companions per BCG is plotted as a function 
of the absolute r' magnitude of the parent galaxy in the left 
panel of Figure [4] This includes all companions down to 20:1 
luminosity ratios. Clearly, this is higher for BCGs than for 
luminosity matched controls. Also, there is a tendency for 
the lower luminosity galaxies (both BCGs and controls) to 
have more 20:1 companions than the brighter galaxies. This 
last, however, may simply be because we probe fainter into 
the luminosity function for the fainter host galaxies. 



3.3 The number and luminosity of companions as 
a function of companion absolute magnitude 

On the left panels of Figure [5] we present the average num- 
ber of companions per parent galaxy (either BCG or con- 
trol) in bins of companion absolute r' magnitude including 
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Figure 5. Left - N c as a function of the companion magnitude for 
BCGs (black) and controls (grey). Right - L c as a function of the 
companion magnitude for BCGs (black) and controls (grey). The 
sample with r p <30 kpc is shown on the top, and the sample with 
r p <50kpc is shown on the bottom. All panels include companions 
down to luminosity ratios of 20:1. 



all companions down to luminosity ratios of 20:1. The his- 
tograms include all companions down to luminosity ratios 
of 20:1 across clusters of all redshifts. The figures clearly 
show that there are few very bright as well as very faint 
companions, with the peak of the distribution of companion 
magnitudes at M r =-20.5. 

The luminosity matched control sample peaks at the 
same companion magnitude as the BCG sample. Although, 
within 50 kpc the number of control companions with 
M r ^-21 is fewer than for the BCGs. N c for companions 
out to only 30 kpc is the same to within errors for BCGs 
and controls; a Kolmogorov-Smirnov test shows no evidence 
that the BCG and control companions were drawn from dif- 
ferent distributions. 

As with N c , the peak in companion luminosity for BCGs 
is higher than for the controls when including companions 
out to 50 kpc. Again, out to only 30 kpc, the two curves are 
the same within errorbars, and a Kolmogorov-Smirnov test 
shows no evidence of a different parent population. These 
results are consistent with the overall N c values reported at 
the beginning of this section. 

Overall, most companion galaxies have moderate mag- 
nitudes of -21. This distribution is most likely a func- 
tion of our parent galaxy sample selection, which peaks at 
M r ^-23, and no BCGs with magnitudes fainter than -21.6 
are in our sample thus no 20:1 companions with magnitudes 
fainter than -18.4 exist in our sample. Since the compan- 
ions are all chosen to have magnitude ratios within 20:1, 
the distribution in companion mass is also expected. Fainter 
companion galaxies certainly exist, but will exist in the less- 
populated bins of the lower-magnitude parent galaxies (see 



Figure [5] right panel). To illustrate this, we calculated N c 
and L c for only the subset of BCGs with magnitudes be- 
tween -22.75 and -23.25. The sample size in this bin is low, 
with only 22 BCGs, however, within the limited magnitude 
range, the number of companions does in fact rise as the 
companion mass falls. Figure [5] reminds us that the bulk of 
BCG companions we explore in this study do not exist as 
dwarf galaxies (at least not down to 20:1), nor as equal-mass 
pairs, but rather in normal L* cluster members. 



3.4 Accounting for projection effects 

Throughout the analysis of the observational data, we have 
studied the number and luminosity of companions which 
are close to the parent galaxy on the sky, not the true 3D 
separation of the companion and parent. This leads to the 
possibility of contamination by foreground and background 
objects, erroneously observed as part of the 2D projection 
of the galaxies on the sky. We can, however, turn to dark 
matter simulations in order to study the importance of these 
projection effects for small separations. We define f3£>, which 
is the fraction of companions which are close in 3D, thus 
allowing for the removal of interlopers. 

Galaxies which are close in three dimensions are those 
likely to merge with the BCG. We attempt to separate the 
potential merging companions from unrelated foreground 



and background galaxies using the Blaizot et al. ( 2005 ) mock 
redshift catalogue, based on the output of the Millennium 
Simulations. After calculating the dynamical friction time- 
scale, we analyse the companions of two samples of simu- 
lated galaxies defined in Section [2 .3| in order to quantify the 
projection effects. Both parent samples match the observed 
data in luminosity, colour, and redshift. They are: (1) All lu- 
minous red galaxies and (2) the CFHTLS BCG-like galaxies 
which are further selected to mimic the cluster core environ- 
ment. 



3.4-1 The dynamical friction time-scale 

There has been much discussion in recent literature about 
the variety of merger time-scales possible. Many authors 
have found that several parameters are important when 
looking at the merging time-scales including gas fractions 



(Lotz et al. 2011) and mass ratios (Williams et al. 2011 
Lotz et al.|2011 ). Based on the same Millennium Simulation 



data, Kitzbichler & White (2008} found that dynamical fric- 



tion time-scales depend on redshift, galaxy mass, as well as 
projected separation. In this work, we do not consider the 
gas mass fraction of the potential mergers, however, we do 
take care to calculate dynamical friction time-scales based 
on a sample of simulated galaxies which has been closely 
matched to our observed data. 



Following Patton et al. (2000), we estimate the merger 



time-scale using the dynamical friction time-scale for our 
samples of high luminosity galaxies. We follow the method 
oflBinney & Tremaine|(| 19871) , assuming circular orbits and 
a dark matter density profile of p(r)oc r -2 . 
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Figure 6. The dynamical friction time-scale: The dynamical friction time-scale plotted as a function of the parent- companion 
separation. The first two panels are for the Millennium simulation data. Left - The real 3D separation and circular velocity is used to 
calculate Tf r i c . Centre - The projected separation is used to calculate Tf r i c . In the Millennium data, the grey points are all luminous 
galaxies and the black dots are the CFHTLS BCG-like galaxies described in the text. Right - Tf r i c calculated for the CFHTLS data, 
using the projected separation and v c = 800km/ s, typical of cluster velocity dispersions. The 126 red dots are companions down to 20:1 
magnitude differences, and the smaller blue dots are the additional 250 (incomplete set of) even fainter companions. 



fric z 



=2. 64x10 s r 2 v c /M ln(A), 



(1) 



where the physical pair separation is r in kpc, v c is the 
circular velocity in kms -1 , M is the mass of the companion 
galaxy in solar masses, and In A is the Coulomb logarithm. 
We use the value of In A ~ 2, measured by Dubinski et al.| 
(1999), since this value agrees with simulations of the orbital 
decay of equal mass mergers. 

We find that even though the galaxies we study likely 
reside at the centre of rich clusters, where the relative veloc- 
ities are often thought to be too high for merging to occur, 
the dynamical friction time-scale is indeed short for a large 
fraction of the companions. This means that close pairs are 
likely to be bound and lead to real mergers. 

In Figure [6] we show that the dynamical friction time- 
scale, Tfric ( |Binney &; Tremaine 1987), is < 0.5 Gyr, for 
close companions within 130 ~ 30 kpc. Ninety percent of 
companions to the luminous galaxies, with projected sepa- 
rations < 30 kpc, have T fric < 0.5 Gyr. All of the CFHTLS 
BCG-like simulated galaxies within 130 ~ 50 kpc have short 
merger time-scales (left panel) . These calculations also show 
that 50% of companions with short merging times are re- 
jected from our measurements of N c and L c as they have 
projected separations > 50 kpc (middle panel). 

For comparison, we calculate Tf r i C for the CFHTLS 
data. We do not have redshift data for the CFHTLS galaxies, 
therefore we use a value of v c =800 km/s, the median value of 
the 3D velocity difference in the Millennium CFHTLS BCG- 
like galaxies and their companions with luminosity ratios 
between 20:1 and 2:1. This value is also consistent with the 
velocity dispersions of massive galaxy clusters, where BCGs 
are located. To find the mass, we turn to the literature for 



characteristic BCG values, van der Marel (1991) determined 



the average M/Lr for elliptical galaxies to be 4.7, and recent 



results from Cappellari et al. (2006) find typical M/Lr of 
5.8 for luminous galaxies with Mr between -22 and -25. We 
take the average of these two values and use a mass to light 
ratio of 5.3. Using this ratio, we convert the luminosity to a 
mass. The trend appears similar to what we found for the 
Millennium data, with a short enough Tf r i C for ^60% of the 
20:1 companions within 30 kpc to merge with their parent 
galaxy in under 0.5 Gyr (Figure [6] right panel). With the 
caveats that the observational data is based on only 126 
BCG companions and our assumption on the velocities, we 
conclude that within the small distances probed here the 
dynamical friction time-scale is short, and most companions 
will likely merge with their parent. 



3.4-2 The fraction of real 3D companions 

Figure [7] shows the fraction of real 3D companions (fsz?) as a 
function of the projected distance. The amount of interloper 
contamination is similarly small for 20:1 and 2:1 companions 
with the closest r p . The simulated data also show that while 
f3£> is only ^0.35 for 20:1 companions at r p ~ 500 kpc, f3D 
jumps to ~0.50 for companions within a projected distance 
of 50 kpc, and to ^0.70 within a projected separation of 
10 kpc. This trend of declining f3£> with projected radius is 
similar for the smaller sample of CFHTLS BCG-like galax- 
ies. For the 2:1 companions, the amount of contamination 
is even less severe, with a greater than 50% purity for pro- 
jected distances lower than 300 kpc. We do not show the 
curve for 2:1 companions to the CFHTLS BCG-like sample 
since the errors are too large due to the small sample size. 
When quoting our final N c in the conclusion, we apply the 
f3D from the statistically most significant 20:1 sample. 
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Figure 7. Interloper contamination: The fraction of real 3D 
companions, far,, is shown as a function of projected separation. 
The amount of interloper contamination is similarly small for 
20:1 (filled grey dots) and 2:1 (open grey dots) companions to 
luminous galaxies as well as for the CFHTLS BCG-like galax- 
ies (filled black dots, includes all 20:1 companions) at small r p . 
fsry falls below 50% within projected separations of 50 kpc for 
20:1 companions, and 300 kpc for 2:1 companions (essentially 
BCG-LRG pairs). Vertical lines highlight distances of 30 kpc and 
50 kpc. 



4 DISCUSSION 

4.1 BCG formation and evolution 

4-1.1 The importance of minor mergers 

Because BCGs are so extreme in size and few in number, 
it is rare to find companions with luminosity ratios of 2:1 
and higher. Our sample of BCGs and controls both show 
a similar value of N c as does the more inclusive sample of 



galaxy surveys studied by Patton & Atfield (2008). Con- 



straining their analysis to a luminosity ratio of 2:1, these 
authors find a value of N c ~0.02 for primary galaxies of mag- 
nitudes -22 < M r < -18. Despite the fact that their sample 
includes less massive galaxies in a range of environments, 
with a range of colours, our study finds slightly larger results 
for BCGs with characteristic magnitudes of M r ~ -22.5, but 
consistent to within errors (N c = 0.04±0.02, including com- 
panions within 30 kpc). 

This analysis adds to the mounting evidence support- 
ing the importance of high-ratio companions (potential mi- 
nor mergers) for galaxy formation. We have shown that the 
luminosity in companions per BCG does not plateau, but 
continues to increase even out to mass ratios of 20:1. There 
is a significant amount of luminosity that exists in the sum 
of lower mass companions - the luminosity in companions 
with luminosity ratios of 4:1 is 1.14±0.28x 10 10 L©, whereas 
it has almost doubled to 2.14±0.31x 1O 1O L0 when including 
20:1 companions. This suggests that minor mergers are at 
least as important as major mergers in building the mass of 
the BCG. This is particularly true for the companions out 
to 50 kpc, with the caveat that they suffer more contamina- 



tion from interlopers. The same is true for the control sam- 
ple. Although the L c in BCG companions rises more steeply 
than for control companions, essentially the luminosity of 
the smaller companions is important for both. Moreover, 
Bluck et al. (2012) have recently shown that for a sample 



of distant luminous galaxies (1.7<z<3) half of the mass re- 
sulting from mergers results from major (1:4) mergers, and 
half from minor (1:100). The building up of the galaxy mass 
through mergers is consistent with the results of Bernardi 
(120091) andlAscaso et al.l (120111), who find BCGs that are 



more local, have larger sizes, but that the BCG luminosity 
profiles are not disturbed, implying a smooth accumulation 
of mass. 



4.1.2 Downsizing 

It is known that more luminous BCGs tend to lie in more 
massive clusters (e.g. [Edge 199l1 |Brough et al.||2008| ), how- 
ever, this does not seem to be dominating the number 
of close companions, since it is the less luminous parent 
galaxies that have the larger number of companions. This 
naturally leads to a fainter integral companion luminosity 
for the brighter parent systems, which is consistent with 
what downsizing predicts. If the light traces mass, and the 
brighter BCGs represent more massive systems, and if the 
larger mass systems formed first, it is consistent that we 
would observe more companions in the smaller systems, 
which could still be actively building up. However, to test 
this hypothesis, the effect of the dwarf galaxies that our 
study may be missing should be considered. First, the lu- 
minosity function is not probed to as faint of magnitudes 
for the brighter BCGs, and second, the larger haloes of 
the brighter galaxies may preferentially cover adjacent faint 
companions. We will explore this in a following paper. 



4.1.3 Evolution with cosmic time 

There are several recent studies regarding merging among 
local samples of large red galaxies, or galaxies specifically 
in groups (e.g. |McIntosh et al.||2008| |Masjedi et al.||2009 



Liu et al. 2009). The most similar sample to the present 
study is that of Liu et ah] ( |2009| , who examined BCGs, 
at redshifts z<0.12 with companions within projected sep- 
arations of 30 kpc, and luminosity ratios of 4:1. They find 
49/515 have close pairs (N c ^0.1), and once folding in mor- 
phological information (tidal tails, etc) deem 18 of these to 
be major mergers (N c ^0.03). To compare our results to 
theirs, we look at N c for BCGs, including only the most lu- 
minous companions, those within 4:1. We find a higher value 
of N c =0.15±0.03 for the moderate redshift sample (^0.1 af- 
ter correcting for interlopers). Although there is a larger 
number of companions in our higher-redshift sample, it is 
important to note that the methodology in both studies is 



very different, and also that Liu et al. (2009) include only the 



red companions. The present study includes all companions 
regardless of colour (though few of the low-ratio companions 
are blue.) However, if the close companions do merge with 
the BCG and are not replaced, one would expect N c to de- 
crease as the universe evolves, which is consistent with the 
difference in N c between our massive cluster galaxies, and 
those of |Liu et al.| pQ09| ) . 
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